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In a Ginzburg-Landau model for parametrically driven 
waves a transition between a state of ordered and one of dis- 
ordered spatio-temporal defect chaos is found. To character- 
ize the two different chaotic states and to get insight into 
the break-down of the order, the trajectories of the defects 
are tracked in detail. Since the defects are always created 
and annihilated in pairs the trajectories form loops in space 
time. The probability distribution functions for the size of the 
loops and the number of defects involved in them undergo a 
transition from exponential decay in the ordered regime to a 
power-law decay in the disordered regime. These power laws 
are also found in a simple lattice model of randomly created 
defect pairs that diffuse and annihilate upon collision. 

Very early in the investigation of pattern-forming sys- 
tems far from thermodynamic equilibrium it has been 
recognized that dislocations in stripe patterns near on- 
set are mathematically closely related to defects in the 
equilibrium xy-raoAel since both systems are described 
by a single complex order parameter A in terms of which 
the dislocations are given by locations of vanishing mag- 
nitude \A\ = 0. One of the fascinating phase-transition 
phenomena studied extensively in the xy-model is the 
Kosterlitz-Thouless transition, which is associated with 
the unbinding of defect pairs |^]. This motivated early 
efforts to identify related phenomena also in pattern- 
forming systems like Rayleigh-Benard convection A 
direct analogy between these systems does not hold, how- 
ever, for a number of reasons. Most significantly, the 
phase transitions occur at finite temperature and are in- 
timately related to the relevance of fluctuations, whereas 
in macroscopic systems like Rayleigh-Benard convection 
the effect of (thermal) noise is negligible in most situa- 
tions H]. Moreover, in the absence of noise there are no 
persistent dynamics in the Ginzburg-Landau with real 
coefficients that describes the equilibrium system. 

More recently, spatio-temporal chaos in pattern- 
forming systems has found considerable attention. A 
number of variations of convection experiments have re- 
vealed various different types of spatio-temporal chaos 
like spiral-defect chaos (e.g. Q]), domain chaos triggered 
by the Kiippers-Lortz instability (e.g. Q), dislocation- 
dominated chaos in binary-mixture convection and 
in electroconvection |Q. Extensive theoretical investi- 
gations have focussed on the complex Ginzburg-Landau 
equation (CGLE) for the complex amplitude of oscilla- 
tions arising in a Hopf bifurcation (e.g. In contrast to 
the Ginzburg-Landau equation describing the xy-model 
the coefficients in the CGLE are complex and make this 



system non- variational. 

Various results have been obtained for the dynam- 
ics of defects in spatio-temporally chaotic systems. In 
binary-mixture convection progress has been made to re- 
construct the full wave pattern from the locations of the 
dislocations ||^ . The probability distribution function for 
the number of dislocations has been measured in electro- 
convection 1^ and found to agree quite well with results 
based on the complex Ginzburg-Landau equation and on 
a simple diffusive model [ p^ . The dynamical relevance 
of dislocations has been demonstrated best so far in work 
that succeeded in extracting the contribution of each dis- 
location to the total fractal dimension of the (extensively) 
chaotic attractor of the CGLE . 

The spatio-temporally chaotic states found in pattern- 
forming systems typically arise from ordered states 
through some transition when a control parameter is 
changed. An interesting question is what actually hap- 
pens when the order of the pattern breaks down. In anal- 
ogy with the melting of two-dimensional crystals [|2 13 
one may expect that defects in the pattern may play an 
important role. In this Letter we present results for a 
transition between two spatio-temporally chaotic states 
in a Ginzburg-Landau model for parametrically excited 
waves [Q. While one state is disordered in space, the 
other retains a stripe-like order despite the chaotic cre- 
ation and annihilation of defect pairs. We characterize 
the break-down of order in terms of the defect dynam- 
ics and find that the transition to the disordered state is 
associated with what one may call an unbinding of the 
pairwise created defects. Tentative results for this un- 
binding transition have been presented earlier in |]l5| , |l4| . 

To illustrate the possible role of the defect dynamics in 
the break-down of order, Fig.|l| presents a space-time dia- 
gram of the defect dynamics in the regime with persistent 
spatial order. The y-location of each defect with positive 
topological charge is shown as a solid circle while the y- 
location of each defect with negative charge is shown as 
a dot. The trajectories of almost all of the defect pairs 
form simple loops in space-time. Thus, while in any given 
system defects are always annihilated in pairs of oppo- 
site charge, in the ordered regime the annihilating defects 
have also been created together. It is in this sense that we 
consider them to be bound pairs. The preservation of the 
stripe-like order in this regime can then be understood 
intuitively, since the dynamics of defects in such simple 
defect loops affect only a very small portion of the sys- 
tem and, moreover, render the system almost unchanged 
after their disappearance. In some cases the space-time 
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loops involve two (see arrow in fig.|]) or possibly three de- 
fect pairs. Then the area of the system that is perturbed 
by the defects is larger, but after their annihilation the 
system is still left in essentially the same state as before. 
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FIG. 1. Space-time diagram of defect trajectories, for 
a = 0.25, c = -l+4j, d = l-|-0.5j, s = 0.2, g = -l-12i, and 
b — 2. Circles and dots denote defects of opposite topological 
charge. 

The orientation and position of the stripe pattern is af- 
fected significantly only between the defects. Therefore, 
one may expect that the destruction of order requires 
that the defects in a pair have to unbind and separate 
arbitrarily far from each other. However, it may also be 
sufhcient if the defects in a given pair are annihilated by 
defects from two other pairs, which in turn are annihi- 
lated by further defect pairs, generating a chain of events 
whose trajectory in space-time is a large loop that spans 
the whole system. Simulations show that in fact even 
in the disordered regime most defects are annihilated by 
their 'own' partner forming a simple loop and most loops 
are small compared to the system size. Thus, if a dis- 
tinction between the two regimes can be made based on 
defect trajectories it must involve more subtle aspects 
like the distribution of defect loops as a function of their 
size. Using a detailed statistical analysis of the defect 
dynamics we show in the following that indeed the sta- 
tistical properties of the defect loops in space-time are 
qualitatively different in the two regimes. 

We consider a Ginzburg-Landau model for paramet- 
rically excited waves in an axially anisotropic system in 
which at the bifurcation point the waves travel only along 
the a;-axis, 

dtA + sd^A = dW'^A + aA + bB + 

c\A\^A + g\B\^A, (1) 

dtB - sd^B = d*\7lB + a* B + bA + 

c*\B\^B + g*\A\'^B. (2) 

Here A and B are the complex amplitudes of right- and 
left-traveling waves. All the coefficients except for the 
group velocity s and the parametric forcing h are com- 



plex. The traveling waves arise in a Hopf bifurcation at 
Or = Re{a) — 0. For < the bifurcation is super- 
critical and the waves exist for > 0. The parametric 
forcing is applied at close to twice the Hopf frequency. 
It therefore induces a resonant interaction between the 
counter-propagating waves at linear order ||l^ , p^ . For 
ttr < and 6^ > -|- of standing waves are excited para- 
metrically that are phase-locked to the forcing. These 
waves correspond qualitatively to the waves excited in 
Faraday experiments. In all of the simulations we use pe- 
riodic boundary conditions in both directions and solve 
(|l|,||) pseudo-spectrally with a fourth-order Runge-Kutta 
scheme using an integrating factor for the linear deriva- 
tive terms. 

Numerical simulations of show two distinct 

regimes of spatio-temporal chaos for the parametrically 
excited waves: a conventional regime in which the spa- 
tial correlation function decays rapidly in an essentially 
isotropic way and a regime of spatio-temporal chaos 
which exhibits a strikingly ordered state in which the cor- 
relation function reveals a strong stripe-like order that is 
also apparent in individual snapshots | p^ , p^ . The transi- 
tion between the two states is discontinuous as indicated 
by a jump in the average number of defects. 

In order to get insight into the role of the defect dy- 
namics in this order-disorder transition we track each de- 
fect from its creation to its annihilation. In particular, 
we identify which defect was annihilated by which de- 
fect. From these data we extract the distribution func- 
tions for the number of defect pairs that are participat- 
ing in a space-time loop and for the spatial and tempo- 
ral extent of the loops. Since the relevant information 
turns out to be in the rare, large loops care has to be 
taken to identify annihilation and creation processes re- 
liably and distinguish them from situations in which a 
defect simply moved relatively far in one time step. In 
our defect-tracking scheme we recursively check the dis- 
tances between defects of equal and of opposite charge. 
If in consecutive time steps two defects of equal charge 
are closer than some threshold value 5^'^ — n ■ Si they 
qualify as a single defect that has moved from one to the 
other position. If more than two defects fall into this 
category the closest defects are taken to be 'continuing' 
defects. Defects that are not continuing defects are candi- 
dates for annihilation and creation events. Among those, 
two defects of opposite charge and closer than a second 
threshold S^'' = n ■ S2 are identified as a pair that was 
annihilated (or created) in this time step. After one step 
of this analysis the same analysis is repeated with new, 
increased values for the thresholds, S'f^'^^'' = {n + 1) ■ Si, 
until all defects have been assigned. 

Fig.|^ shows the relative frequency of loops consisting 
of at least n defect pairs. The results in figs.§,||,H are 
based on 8000 timesteps {dt = 0.5) with an average num- 
ber of 7000 defects at any given time in the disordered 
regime. In the ordered regime (b > 0.7) very few loops 
contain more than 5 defect pairs and the distribution de- 
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cays essentially exponentially. In the disordered regime 
(b < 0.625), however, the number of loops with many de- 
fects is greatly increased and the decay of the distribution 
function is only algebraic with an exponent of a « 1.5. 




Number n of Defects in Loop 

FIG. 2. Relative frequency of loops made up of at least n 
defect pairs. Parameters as in fig.^ except for b. System size 
L = 1088 in the disordered regime; L — 272 and L = 136 for 
b — 0.7 and b = 1.0, respectively. 

Since the defect motion essentially affects only the 
(vaguely defined) part of the system between the defects 
a better indicator for the expected loss of order are the 
spatial extents Aa; = x^ax ~ ^min and Ay = ymax ~ Umin 
of the defect loops in the x- and ^/-direction, respectively. 
Here Xmin,max and ymm^max rcfcr to the minimal and 
maximal values of x and y in any given loop. Note 
that these values need not be obtained at the same time. 
Thus, in principle a small loop could still yield large Ax 
or Ay if it traveled. 
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FIG. 3. Relative frequency of loops with spatial extent in 
the a;-direction at least Aa; . Parameters as in fig.pl 



Similar to fig.^j, fig-H shows the relative frequency of 
loops with size in the a;-direction larger than Ax. In the 
ordered regime the decay is again very rapid and there 
are essentially no loops with Ax larger than 10, which 



is of the order of one wavelength. This may indicate 
some pinning of the defects by the pattern in that it may 
restrict their motion to a predominantly climbing mo- 
tion. In the disordered pattern, on the other hand, a 
noticeable number of loops is of the order of the system 
size L — 1088 and even larger (i.e. wrapping around 
the system due to the periodic boundary conditions). In 
a simple interpretation these events would be associated 
with a persistent change in the average wavevector of the 
pattern. Again, the distribution functions are exponen- 
tial in the ordered regime and exhibit a power law in the 
disordered regime with exponent /? w 3. The distinction 
between the regimes is not quite as striking in the spatial 
extent Ay in the y-direction. Even in the ordered regime 
the loops can reach a size of Ay « 100, while the loops in 
the disordered regime extend to sizes of Ay sa 1000. The 
distribution is still exponential in the ordered regime and 
appears to be power-law in the disordered regime. How- 
ever, the measured exponent increases from about 2.8 to 
4 as 5 is increased from 0.4 to 0.625. 




FIG. 4. Relative frequency of loops with lifetime at least 
At. Parameters as in fig.^. 

The duration At of the loops is also of interest and 
the corresponding relative frequencies are shown in fig.^. 
The exponential and power-law character of the distri- 
butions are quite clear in the respective regimes with the 
exponent of the power-law 7 « 2.7. 

To give further support for the existence of power laws 
in the distribution functions in the disordered regime and 
to get more insight into their origin we have investigated 
a simple lattice model of the defect dynamics, which is 
based on the observation that the single-defect statis- 
tics in the disordered regime show the same signature as 
those of the defect chaos in the CGLE. This type of dis- 
tribution has been shown to arise if the defects behave as 
random walkers that are annihilated upon colliding with 
any other defect of opposite charge ||l^ . 

The results of our implementation of the simple lattice 
model are shown in fig.^ With a probability p, ran- 
dom walkers of opposite charge are created pairwise at 
the same randomly chosen site of a square lattice. They 
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interact only with walkers of opposite charge and an- 
nihilate upon contact. Fig.|| shows the loop statistics 
for these walkers in a system of size L — 1600 and a 
probability of creation p = 0.000016 (thick lines) and 
p = 0.0001 (thin lines). All measured quantities, i.e. 
number of defects in a loop and the spatial as well as 
the temporal extent of the loops, show power-law behav- 
ior for large loops. The simplicity of the lattice model 
suggests that the power laws are expected to arise in a 
much wider class of spatio-temporally chaotic systems 
including the CGLE. The exponents are measured to be 
a = 1.6, (3 = 2.9, and 7 = 2.4, respectively. Thus, even 
the values of exponents agree quite well with those ob- 
tained in the simulations of (|l|,||). 
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FIG. 5. Simulations of the lattice model. Thick lines de- 
note p = 0.000016, thin hues p = 0.0001. L = 1600. 

In conclusion, using numerical simulations of two cou- 
pled complex Ginzburg-Landau equations we have in- 
vestigated the break-down of order in a transition from 
a non-equilibrium, chaotic stripe phase to a disordered 
phase. To get insight into the way the order breaks 
down we have analyzed the trajectories of the disloca- 
tions in the pattern. In particular, we have determined 
the statistics of the loops formed in space-time by chains 
of creation and annihilation events of oppositely charged 
defect pairs. The order-disorder transition is related to 
a significant increase in the number of loops that extend 
over the whole system, which we associate with the un- 
binding of defect pairs. More precisely, the decay of the 
loop distribution function changes from exponential to 
algebraic in that transition. The algebraic decay is also 
found in a simple lattice model of diffusing and anni- 
hilating defects, with the exponents agreeing quite well 
with those found in the Ginzburg-Landau equations. The 
agreement between the Ginzburg-Landau equations and 
the lattice model suggests that the power laws may be 
found in a wider class of defect-chaotic systems and in 
particular also in the single complex Ginzburg-Landau 
equation. 

The single-defect statistics obtained in the disordered 
state and in the lattice model have also been found in 



experiments in electroconvection and in thermal con- 
vection in an inclined layer |l8| . It would be exciting to 
study also the loop distribution functions in these sys- 
tems. The amount of detailed data necessary to deter- 
mine the loop statistics represents presumably a chal- 
lenging task in the current set-ups. 

The simplicity of the lattice model suggests that the 
power laws may also be amenable to an analytic ap- 
proach. A further interesting question is whether the 
lattice model can be extended to obtain also a transition 
from power-law to exponential decay. Introducing an at- 
tractive interaction between the defects would seem to be 
a natural choice. It is not clear whether the interaction 
should be short- or long-range. 
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